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Abstract. We modify the conventional curvature radiation (inverse Compton scatter- 
ing) + synchrotron radiation polar cap cascade model by including the inverse Compton 
scattering of the higher generation pairs. Within the framework of the space-charge- 
limited-flow acceleration model with frame-dragging proposed by Harding & Muslimov 
(1998), such a full polar cap cascade scenario can well reproduce the L 7 oc (L^) 1 ^ 2 
and the L x ~ 10 _3 L s( j dependences observed from the known spin-powered pulsars. 
According to this model, the "pulsed" soft ROSAT-band X-rays from most of the mil- 
lisecond pulsars might be of thermal origin, if there are no strong multipole magnetic 
components near their surfaces. 

INTRODUCTION 

Eight and 35 spin-powered pulsars have been also detected in 7-ray and X- 
ray bands, respectively. Despite their great diversity of emission features, the 
luminosities of these pulsars seem to obey the empirical laws L 7 oc (L s d)^ 2 oc B/P 2 
(Thompson et al. 1997), and L x ~ 10" 3 L sd oc B 2 /P 4 (Becker & Triimper 1997), 
where L s d is the spin-down luminosity of the pulsar. The spectra of the 7-ray 
emission and the X-ray emission from most of the pulsars are non-thermal, while 
full surface thermal emission components are identified from Vela, Geminga, PSR 
1055-52 and PSR 0656+14, and possible hot polar cap thermal emission components 
are reported from PSR 1929+10 and PSR J0437-4715. 

Two competing models for pulsar high energy emission, i.e., the polar cap mod- 
els (Daugherty & Harding 1996; Sturner, Dermer & Michel 1995) and the outer 
gap models (Cheng, Ho & Ruderman 1986; Zhang & Cheng 1997) were proposed. 
Canonical polar cap cascade models involve the curvature radiation (CR) or inverse 
Compton scattering (ICS) of the primary particles and the synchrotron radiation 
(SR) of higher generation pairs. Here we modify such a cascade picture by including 
the ICS of the higher generation pairs, which is important since it usually occurs in 
the resonant regime. A more detailed presentation of this study is shown in Zhang 
& Harding (1999). 



THE MODEL 



The "full-cascade" picture 

The "full-cascade" scenario is: primary particles accelerated from the inner gap 
emit primary 7-rays via CR or ICS, these 7-rays will pair produce in strong mag- 
netic fields. The secondary pairs have non-zero pitch angle with respect to the field 
lines. The perpendicular energy of the pairs will go to high energy radiation via SR, 
and the parallel energy of the pairs will also convert to radiation via ICS with the 
soft thermal photons from either the full neutron star surface or the hot polar cap. 
Under the condition of 7 > 7 rcs = 48-Bi 2 Max[(l — /3/i S)inax )T Sj6 , (1 — Pfi>h 
(i.e. the "resonant scattering condition", Dermer 1990), the efficiency of converting 
particles' kinetic energy to radiation by ICS is almost 100%, so that the total high 
energy emission luminosity is approximately the polar cap particle luminosity. Here 
T is the temperature of the soft photons, and /i max is the cosine of the maximum 
scattering angle. Note that two components, i.e., a soft full surface thermal emis- 
sion (denoted by 's') and a hard hot polar cap thermal emission (denoted by 'h'), 
are adopted. 

The basic ingredients in constructing an analytic description of such a full cascade 
scenario are (1) the energy distribution in SR (perpendicular) and ICS (parallel) 
branches, and (2) the recursion relations between different generations. For the 
former, r)\\ = 7i,||/7i = [1 + (7 2 + i — 1) sin 2 ^ kB ] -1 ^ 2 and r/j_ = 1 — r)\\ are the energy 
portions for the parallel (ICS) and the perpendicular (SR) branches, respectively, 
where % and 7^1 are the total and parallel Lorentz factors of the i-th generation 
pairs, and 9 kB is the impact angle between the photon and the magnetic field line. 
With this, one can get the reduction factor of the typical photon energies for the 
adjacent generations, e.g., k S r = e i+liSR /ei = (3/4)x, and k IC s = ei+i.ics/e* = r)\\B' e , 
where B' e = B e /B cri , and x is the parameter to describe the 7 — B pair production 
process. The photon escaping energy is -E 7)CSC (nPSR) ~ 2.0GeVB~l 2 P 1 ^ 2 r^ 2 x 1/16 
for normal pulsars, and E 7i e SC (msPSR) ~ 73GeV B~l P^ 2 r~£ 2 x 1/12 for millisecond 
pulsars, where X 1/16 = x/(l/16), and X 1/12 = an d r e $ is the emission 

height in units of 10 6 cm. Given the typical primary photon energy E (which is 
model-dependent), we can then get some non-integer generation order parameters 
(Lu et al. 1994; Wei et al. 1997), e.g. Csr = l2 ^§^p + 1 ( for P ure SR generations), 

and (jcs = ^^^7^ + -'- (f° r P ure ICS generations), which can describe the complex 
cascade process analytically. 

Harding & Muslimov acceleration model 

Harding & Muslimov (1998) has improved the space-charge-limited flow accel- 
eration model (Arons & Scharlemann 1979) by incorporating the frame-dragging 
Eu , upper and lower pair formation front and both the CR and ICS of the primary 



electrons. It was found that a stable accelerator is located at an effective "radius" 
of Re ~ (1.5 — 2)R for normal pulsars when ICS energy loss is less than that of 
CR, since the ICS of the upward versus downward particles with the soft ther- 
mal photons are anisotropic due to different geometries. The typical length of the 
CR-controlled gap is S c ~ 4.8 x 10 4 cm£?~ 4 2 7 ' P^ 7 Res (cos a)" 3 / 7 , and the typical 

Lorentz factor of the primary particles is 70 = 4.7 x 10 7 B~^ P 1 ' '"' 'r^fi (cos a) 1 ' 17 , 
where a denotes the inclination angle of the neutron star. Thus the typical energy 
of the primary photons is E ~ 33.2(GeV)S~i2 7 P~ 1//14 re 7 6 //14 (cosQ;) 3 / 7 , with which 
one can get explicit expressions for the generation parameters, and complete the 
analytic description of the full-cascade model. 

LUMINOSITY PREDICTIONS 

Gamma-ray luminosity 

An interesting feature of the Harding & Muslimov model is that 70 is insensitive 
to pulsar parameters, so that the polar cap luminosity L pc = r y mc 2 N p is roughly 
proportional to (L^) 1 / 2 . One advantage of the full cascade model is that the ICS 
branches can convert the "lost" parallel kinetic energies of the particles also to 
radiation, so that the total high energy luminosity (mainly 7-ray luminosity) is 
also roughly proportional to (Lgd) 1 / 2 . More specifically, the model predicts 

L 7 (full) ~ L pc = 5.4 x 10 31 erg • s" 1 B^p- 13 / 7 ^ (cosa) 8/7 (la) 
= 1.7 x 10 16 -B-g 'P^ VJ? '(cos afr< '(L sd )^ 12 . (lb) 

This nearly reproduces the observed L 1 oc (Lgd) 1 ^ 2 feature. 

Thermal X-ray luminosity 

There are two thermal components in a pulsar's X-ray spectrum (though they 
might be buried under the non-thermal component). For the full surface thermal 
component, we have adopted a simple rough "standard" cooling model, which is 
not inconsistent with the observations. 

For the hot polar cap thermal component, we treated it with the self-consistent 
polar cap heating in the Harding & Muslimov model. Since the flow is space-charge- 
limited so that the deviation of local charge density (p) from the Goldreich- Julian 
density (pgj) is small, the backflow particle luminosity should be only a small 
portion (a factor of | p — pgj I /2pcj) of the polar cap luminosity, which reads 

L x , th , max ~ 5.9 x 10 29 erg ■ s" 1 B^P^^n ' (cos a) 5/7 > (2) 

and the maximum polar cap temperature (assuming an area of nr 2 , where r p = 
9 pc R = 1.45 x 10 4 P- 1 / 2 cm) is 
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FIGURE 1. Observation versus theory: (a) 7-ray luminosities (b) X-ray luminosities. 



T pc , max = 2.0 x lO^S^P-^r-J^Coosa) 5 / 28 . 



(3) 



Non-thermal X-ray luminosity 

The SR spectra of all SR branches can not get down to X-ray energies as those 
observed by ROSAT and ASCA, since there is a low energy cut-off at the blueshifted 
local resonant frequency (corresponding to the transition between the ground and 
the first excited Landau levels). Another advantage of the full cascade scenario is 
that, the ICS branches, which were neglected in the canonical cascade model, can 
naturally give a non-thermal component extend to the X-ray band. For normal 
pulsars, the maximum non-thermal luminosity below a certain energy E c could be 
estimated as 
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where Cics,k — '"fo^tic^ + 1 is the number of pure ICS generations for the typical 
energy of the fc-th SR generation, E k = E k,^ r , to reduce to the escaping energy 
E csc , Vc,k,j = Icllkj (but = 1 when 7 C > 7 fej -), 7 C = E c /[(1 - f3/j,)2.8kT], and 7^ = 
(e /2)/tg^ 1 /t^g?7||. For millisecond pulsars, a slightly different formula is adopted 
(see Zhang & Harding 1999). When calculating the X-ray luminosity within a 
certain band, E cl < E < E c2 , we then have L AE = L xMh (E c2 ) - L x , nth (E cl ). 



RESULTS AND DISCUSSIONS 



The observation versus model prediction of the 7-ray and X-ray pulsars are shown 
in Fig.l and Fig.2. For the X-ray luminosities, three components (the non-thermal, 
the full surface thermal and the hot polar cap thermal components) are taken into 
account. 

An obvious conclusion is that the full polar cap cascade model within the 
framework of Harding & Muslimov acceleration model can both reproduce the 
L 7 oc (Lgd) 1 / 2 and L x (ROSAT) ~ lCT 3 L s d feature simultaneusly. The former was 
not done by the outer gap model which interprets non-thermal X-ray emission as 
the SR of the downward cascade from the outer gap (Cheng & Zhang 1999; Zhang 
& Cheng 1997). In our model, we also compare the non-thermal X-ray luminosity 
with the luminosities of the two thermal components. It is found that for middle- 
aged pulsars such as Vela, Geminga, PSR 1055-52, and PSR 0656+14, the full 
surface thermal luminosity is comparable to the non-thermal one, so that such a 
component should be detectable from these pulsars. Such a feature is actually 
observed. For the hot polar cap thermal component, our model shows that it is de- 
tectable in relative old pulsars such as PSR 0950+08 and PSR 1929+10, although 
the non-thermal component is also detectable. The outer gap model actually only 
predicts pure thermal components in these pulsars, since the thick outer gap does 
not exist in their magnetospheres (Cheng & Zhang 1999). For the millisecond pul- 
sars, our model predicts that the thermal emission from polar cap heating is the 
dominant component in the "pulsed" ROSAT-band spectra, while the outer gap 
model predicts non-thermal emission, though they assume strong multipole mag- 
netic fields near the surfaces of the millisecond pulsars. Future observations and 
spectral analyses can distinguish the riveling models. 
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